Acoustic Emission and Shape Memory Effect in the Martensitic Transformation 
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Acoustic emission signals are known to exhibit a high degree of reproducibility in time and show 
correlations with the growth and shrinkage of martensite domains when athermal martensites are 
subjected to repeated thermal cycling in a restricted temperature range. We show that a recently 
introduced two dimensional model for the martensitic transformation mimics these features. We 
also show that these features are related to the shape-memory effect where near full reversal of 
morphological features are seen under these thermal cycling conditions. 

PACS numbers: 81.30.Kf, 75.60.Ej, 45.70.Ht 



Acoustic emission(AE) is an important nondestructive 
technique that is sensitive to the microstructural changes 
taking place inside the system. It is widely employed in 
the detection of earthquakes of small magnitudes , un- 
derstanding and mapping nucleation events of fracture 
in seismology Q , and precursor effect of fracture Q • Its 
nondestructive character finds many applications in in- 
dustry as well Q. Acoustic emission is attributed to the 
sudden release of stored strain energy, although the de- 
tails of the mechanism are generally system specific. For 
instance, the difference between the nature of AE signals 
in thermal and athermal transformations has been used 
to characterize the nature of a first order transformation 
In the case of athermal martensites, the kinetics of 
the transformation is not controlled by thermal fluctua- 
tions which implies that the system is caught in a local 
minimum. It jumps over the barrier only when subjected 
to changes in external parameters such as temperature or 
stress |1 H HUH El- Thc AE emitted during marten- 
sitic transformation often shows unusual and interesting 
statistical properties, one that has attracted considerable 
attention is the power law statistics of the AE signals 

f, 0| , as it is reminiscent of self-organized criticality 
Much less known, but perhaps equally intriguing is 
that AE signals can show a high degree of statistical cor- 
relations in time (reproducibility) and its correlations to 
the growth and shrinkage of microplates 0, 0, ^| . These 
results based on careful experiments on single crystals 
of CuAlZn alloys under repeated thermal cycling in a 
restricted domain of temperatures 0, IT^ | are well sum- 
marized in Figs. 10-12 of Ref. To the best of our 
knowledge, these results have not been explained so far. 
Explaining these features is also crucial to the under- 
standing of the shape-memory effect of martensites . 
Moreover, as far as we are aware, there are no models 
of shape-memory effect which demonstrate the reversal 
of morphology of martensite domains under thermal cy- 
cling. The purpose of this paper is to devise a model 
which explains these correlations and provide a model 
for the shape-memory effect. 

Martensitic transformation is considered as an atypical 
first order transformation as it exhibits some features of 
second order transition such as precursor effect and 



the scale free power law distribution of AE signals, again 
a signature of critical fluctuations. On cooling from a 
higher symmetry parent phase, nucleation of thin plate 
like product domains with internal twinned structure is 
seen. Nucleation is known to be athermal usually oc- 
curring at defects [16j . Long range strain fields resulting 
from the difference in the unit cell structure between the 
parent and product phases tend to block the transfor- 
mation. The athermal nature of the transformation also 
means that the amount of transformed phase is fixed for 
a fixed quench and further undercooling is required for 
an additional growth. For the same reason, the transfor- 
mation occurs in a broad range of temperatures and it 
exhibits hysteresis under thermal cycling. 

The power law statistics of the avalanches has been 
modeled using random field Ising models [17j | where the 
extent of quenched disorder plays a crucial role. How- 
ever, as the power law is seen only after repeated thermal 
cycling, Vives et al 0,0] conclude that the transforma- 
tion induced disorder, rather than quenched disorder ( 
defects ) is responsible for these avalanches. Recently, 
a model which generates dynamical disorder during the 
transformation has been shown to reproduce the power 
law statistics ^lj- We show that this simple 2D phe- 
nomenological model nas the right ingredients to ex- 
plain the correlated AE signals and the shape-memory 
effect 0. 

The model uses the deviatoric strain as the principal 
order parameter. The effect of the bulk and shear strains 
are included only in a phenomenologically way through a 
long range interaction term (see below). In terms of the 
two components of the displacement field, u x and u y , the 
deviatoric strain is given by e = (^r — ^-)/\/2. 

The scaled free-energy, written entirely in terms of e(r), 
is the sum of local a free energy Fl and an effective long- 
range term, Fi r , that describes transformation induced 
strain-strain interaction. Fl is given by 
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where both a and a are in scaled form. Here, /j(e(f)) 
is taken as the usual Landau polynomial for a first order 
transition given by fi(e(r)) = §e(r) 2 — e(r) 4 + ^e(r) 6 , 
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where r = (T — T C )/(T — T c ) is the scaled temperature, 
T is the first-order transition temperature, and T c is 
the temperature below which there are only two degen- 
erate global minima e = ±e m . We mimic the nucleation 
of the martensite domains occurring at localized defects 
sites by an inhomogeneous stress field, a(r) Recent 
studies have shown that an effective long-range interac- 
tion between the transformed domains results due to the 
coupling with the other components of the strain order 
parameter |l5l Il9| . However, as our approach is phe- 
nomenological, we attempt to guess a simple form that 
respects the invariance of /{(e) under e — > — e. Follow- 
ing Wang and Kachaturyan 20] (see their Eqn (11)), one 
simple form that satisfies this requirement when repre- 
sented in the Fourier space can be written as 



(2) 



F lr {e} = / dkB(k/k){e 2 (r)} k {e 2 (r)} k , 



where {e 2 (^)}fe and {e 2 (^)} fe « is the Fourier transform 
of e 2 (r) and its complex conjugate respectively, and 
B(k/k) is an appropriate kernel. Further, noting that 
the favourable directions of growth of the twins are along 
[11] and [11] directions, we find that the functional de- 
pendence of the kernel on k x and k y takes on a rela- 
tively simple form. We also use the fact that [10] and 
[01] are unfavourable directions for the growth. A sim- 
ple choice of B(k) having these features is B(k/k) = 
— \(3Q{k — A)k^ky, where k x and k y are the unit vectors 
in the x and y directions, and (3 is the strength of inter- 
action. The step function 9(k — A) has been introduced 
to impose a cutoff on the range of the interaction. (See 
Ref. 21].) Equation 2 also has the property that as 
the transformation proceeds, the growth of the domains 
transverse to [11] and [11] directions are hindered. In- 
deed, the geometrical picture of the kernel in real space 
is similar to Kartha et al 0] . 

Acoustic emission during the transformation suggests 
that inertial effects are important j2^. This is included 
through the kinetic energy of the displacement fields 
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where p is the mass density. Experiments show that 
the AE activity depends on the acceleration of the mi- 
croplates |l3|. Since, the parent-product interface moves 
in the parent phase, it is associated with dissipation 
which we represent by the Rayleigh dissipative functional 
|22| . Again, we represent the dissipative functional en- 
tirely in terms of the deviatoric strains |23[ . 



(4) 



This is consistent with the fact that shear and bulk 
strains are known to equilibrate rapidly and hence do not 
contribute. The equations of motion for the displacement 



fields are calculated using the Lagrangian L = T — F 
where F is the total free energy [23| through 
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After scaling out p and a, the equation of motion for the 
deviatoric strain e(r, t) is given by 



dt 



+ 2e{r,t) / dkB(k/k){e 2 (k,t)} k e 



ik.r 
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where (3 and 7 now refer to the rescaled parameters. 
There are three adjustable parameters in the model, 
namely, the scaled temperature r, the strength of the 
long-range interaction (3 and that of dissipation 7. The 
effect of increasing 7 is to reduce the twin width and that 
of [3 is to reduce the width of the martensite domains 
transverse to the habit plane directions. 

Simulations are carried out by discretizing Eq.(6) on 
a, N x N grid ( ie., Ax = 1) using the Euler's scheme 
with periodic boundary conditions. The time step used 
is At = 0.002. ( Actually, we see convergence of solu- 
tions even for At = 0.005.) A pseudo-spectral technique 
is employed to compute the long-range term [Til ] . Com- 
putations are carried out for several values of parameters. 

Our earlier study has shown that the model exhibits 
hysteresis under thermal cycling during which the en- 
ergy dissipated R(t), (— ~^dE/dt) occurs in the form 
of bursts 0. Using Eq.(4), we compute R(t) during 
repeated thermal cycling in a small temperature inter- 
val as in experiments. We represent the stress field 
due to the random distribution of defects by cr(r) = 
Yjj max cro(fj)exp[— \f— rJ| 2 /C 2 ]i where fj refers to the co- 
ordinates of the randomly chosen defect sites and j max 
their total number. (Jo(rj) is taken to be uniformly dis- 
tributed in the interval [-0.3,0.3]. Using the above ini- 
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FIG. 1: (a) Area fraction <f> during thermal cycling in the 
interval r =-5 to 4 for cycles 1 (*) , 2 (<) and 3 (+) . (b) 
Hysteresis for the full cycle from the Austenite to Martensite 
phase and back. (•) is the starting point of the cycles. 

tial configuration, we cool the system from the austen- 
ite phase. Then, at an appropriate point ( shown as • 
in Fig. 3b) in the full cycle, we subject the system to 
repeated thermal cycling in a small temperature range 
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Tmin = — 5 and T max — 4. The starting configuration of 
the martensite domains for the small cycle is chosen in 
a way that, it contains a fair number of martensite do- 
mains, typically, when the area fraction 4> ss 0.6. For the 
cooling cycle, the final configuration during heating cy- 
cle is taken as the initial configuration. Calculations are 
performed for a range of parameter values of 50 < (3 < 10 
and 5 < 7 < 1. As the general features remain the same, 
we report here the results for f3 = 25 and 7 = 4. Other 
parameter values are A = 0.2, £ = 1 and N — 256. Dur- 
ing the first few cycles, the loops drift slightly in the 0— r 
plane which mimics the training period known in exper- 
iments, eventually stabilizing after several cycles (Fig. 
la). (In experiments also, the training period is typically 
the same for CuAlZn @, but could be much larger in 
some other alloys.) A typical plot of <f> vs. t for the first 
three cycles is shown in Fig. la. 

The energy dissipated R(t) stabilizes after the training 
period. A plot of R{t) for several forward and reverse cy- 
cles ( fifth to eight) after stabilization is shown in Fig. 2. 
It is clear that the energy bursts (which mimic the AE sig- 
nals), as in experiments |slll3||. exhibit a near repetitive 
pattern in time (temperature) during successive heating 
and cooling parts of the cycles. 
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FIG. 2: Repetitive nature of R(r) for cycles 5 to 8. 



In order to establish a correspondence between R(t) 
and the changes in the spatial configuration of the 
martensite domains, we have simultaneously monitored 
the morphology over several cycles. During the first few 
cycles, the morphology changes considerably even as the 
system returns to nearly the same point in the <fi — r dia- 
gram at the end of each cycle. Most changes occur during 
the first three or four cycles as can be seen from Figs 3a 
and 3b. Figure 3a is the starting morphology for the first 
cycle, while Fig. 3b that obtained at the end of four cy- 
cles. As can be seen, all the curved twin interfaces have 
been rendered straight after going through four cycles. 
The snapshots of the morphology during one such stabi- 
lized cooling and heating cycle, the fifth one, at selected 
intervals is displayed in Fig. 4. It is clear that during 
heating, the martensite plates shrink and some even dis- 
appear. However, during cooling these domains reappear 
and the eventual morphology at the end of the cycle is 
practically recovered on returning to the starting point 
on the (<p, t) diagram. As can be seen, the final configu- 
ration obtained during the fifth cycle, Fig. J^d can be seen 
to be practically the same as Fig. 3b which is the initial 



FIG. 3: Morphology of the initial configurations for (a) Cycle 
1 (b) Cycle 4. Grey cells correspond to the austenite phase, 
and black and white to the two martensite variants 

configuration for the fifth cycle. Under further cycling, 
the end configurations change very little, consistent with 
that observed in experiments [|J . 



t=-2.0 (a) x=-5.0 (d) 




FIG. 4: Sequential morphological snapshots for r = -2.0, 4.0, 
-0.5 and -5.0 during the fifth cycle. The initial configuration 
for the cycle is shown in Fig. 3(b). 

Having established the repetitive nature of the energy 
bursts and the reversal of the morphology under thermal 
cycling, we need to explain the role of training cycles 
in inducing this correlated behavior shown in Figs. 2 
and 4. We have argued earlier that the model gener- 
ates a huge number of transformation induced metastable 
states due to the fact that we have including threshold 
dynamics, slow driving ( the rate of cooling or heating), 
dissipation and a fast relaxation mechanism. The spiky 
nature of R(t) which mimics the bursts of acoustic emis- 
sion was attributed to the following. As the temperature 
changes, the increase in the driving force arising from the 
decrease in temperature , goes in surmounting these local 
thresholds and another part goes in to the growth of the 
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platelets, and the rest is dissipated in the form of bursts 
of energy (due to advancing interface). First, we note 
that the initial state for the repetitive small thermal cy- 
cles, Fig. 3a, is a configuration reached by cooling from 
the austenite phase. Further, a comparison of Fig. 3a 
with the final morphology reached after four small ther- 
mal cycles, Fig. 3b, reveals that the initial state Fig. 
3a is a relatively shallow metastable state. This can be 
seen from the fact that the sizes of the martensitic do- 
mains of Fig. 3a, are smaller than that in Fig. 3b (or 
Fig. 4d) and the twin interfaces of Fig. 3a are rough and 
considerably curved compared to Fig. 3b. (Note that 
the austenite-martensite interfaces are nearly straight). 
Such configurations are generally expected to have higher 
energy compared to straighter ones. The repetitive pat- 
tern of the energy bursts during successive cycles after 
the training period is an indication that the system tra- 
verses through the same set of metastable states. Dur- 
ing the first few cycles, the free energy landscape is so 
modified that it smoothens out the high energy barriers 
corresponding to the rough curved twin interfaces in Fig. 
3a with very little change in the area fraction. A cru- 
cial role in smoothening process is actually played by the 
long-range interaction term, as the growth (shrinkage) of 
a martensite domain is influenced by the configuration 
of rest of the domains. Computation of the free energy 
Fi r arising from the long-range interaction between the 
domains shows that it actually becomes more negative, 
saturating after the first few cycles. This leads to a re- 
duction in the local free energy, Fl, as well. The net 
effect is to create a deeper set of metastable states for 
the system to circulate for the later cycles. Within one 
such stabilized cycles, say the fifth, the starting config- 



uration (Fig. 3b) has the lowest free energy reaching a 
maximum at the end of a heating cycle (Fig. 4b) . There- 
after it decreases during cooling. The fact that after first 
few cycles, the system traverses through a deep set of 
free energy configurations is actually reflected in Fig. 4. 
Note that the morphology at t = 4.0 is similar to that at 
t = —0.5 except that the widths have increased at the ex- 
pense of the austenite regions consistent with the increase 
in the well depth stated above. The growth in the width 
of martensitic domains as we decrease the temperature 
is surprisingly similar to that observed in experiments. 
(See Fig. 12ofRef.0.) 

In conclusion, our model exhibits the repetitive bursts 
of energy under successive thermal cycles in a small tem- 
perature interval as observed in experiments on acous- 
tic emission. More importantly, these bursts of energy 
are shown to be correlated to the growth and shrinkage 
of martensite plates. The near full reversal of the mor- 
phology during successive cycles after the training period 
can be traced to the influence of the long-range interac- 
tion term. During the training period the long-range 
term has a tendency to smoothen out higher energy bar- 
riers in the free energy landscape. This in turn induces a 
transformation pathway along a unique set of low energy 
metastable configurations. As far as we know, this is the 
first time that the influence of training cycles on shape- 
memory effect has been elucidated. We expect that our 
analysis provides a good insight into the shape memory 
effect which finds immense applications in a variety of 
areas from mechanical actuators to bio-medical applica- 
tions 0- To the best of authors knowledge, this is the 
first model which shows near full reversal of morphology 
under thermal cycling. 
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